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ABSTRACT 

We investigate the contribution made by active galactic nuclei (AGN) to the high- 
redshift, luminous, submillimetre (submm) source population using deep 2 
mJy/beam) Large Apex Bolometer Camera (LABOCA) 870/im observations within 
the William Herschel Deep Field. This submm data complements previously obtained 
Chandra X-ray data of the field, from which AGN have been identified with the aid 
of follow-up optical spectra. From the LABOCA data, we detect 11 submm sources 
(based on a detection threshold of 3.2ct) with estimated fiuxes of > 3 mJy/beam. Of 
the 11 identified submm sources, we find that 2 coincide with observed AGN and that, 
based on their hardness ratios, both of these AGN appear to be heavily obscured. We 
perform a stacking of the submm data around the AGN, which we group by estimated 
Nh column density, and find that only the obscured {Nh > 10^^ cm~^) AGN show 
significant associated submm emission. These observations support the previous find- 
ings of Page et al. and Hill & Shanks that obscured AGN preferentially show submm 
emission. Hill & Shanks have argued that, in this case, the contribution to the ob- 
served submm emission (and thus the submm background) from AGN heating of the 
dust in these sources may be higher than previously thought. 

Key words: galaxies: high-redshift quasars: general submillimetre: galaxies X-rays: 
galaxies. 



1 INTRODUCTION 

It is thought that the production of submm emission in the 
. luminous galaxy population observed at ~ 850 ^m is pre- 
' dominantly driven by star-formation (e.g . I S mail et al.ll2004l : 
I [Alexander et al.ll2005al : iPope et al.ll20od '). However, there is 
■ increasing observational evidence for a significant contribu- 
' tion from AGN to the heating of dust, and hence the produc- 
, tion of submm emission, in these galaxies. Certainly some 
fraction of submm sources are observed to emit hard X-rays 
and there is still the possibility that the observed X-rays are 
reflected or scattered, which would mean they could have 
a significantly higher luminosity at 10-30 keV energies than 
expected. 

Obscured AGN (i.e. with Nh > 10^^ cm"^) in partic- 
ular appear to show high submm fluxes, as identified both 
from statistical techniques (|Hill fc Shanks|[2011bl : [Lutz et al.l 

I2OI0I ) and targeted submm observations of known QSO s 
(e.g. ICoppin et all l2008bl : iMartmez-Sansigre et all [200i). 

iPage et al.1 ( 20041 ') presented samples of obscured and un- 
obscured QSOs and suggested that there may be a differ- 
ence between these two classes in the submm, with the 
obscured sources showing higher submm fluxes (see also 
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iMainieri et al. I I2OO5I for a fainter example of Type 2 QSO). 
This is a crucial observation since it implies that the unifled 
AGN model - where X-ray obscured QSOs are interpreted 
as being viewed at a different angle than unobscured QSOs - 
may be incomplete. Instead, X-ray obscured QSOs may rep- 
resent an earlier stage in QSO evolution, where the black- 
ho le formed w i thin a dusty starburst galaxy, as suggested 
bv lPage et al.l (|2004l V This would account for an instrinsic 
difference between the submm properties of obscured and 
unobscured AGN. 

Recently, ICarrera et al.1 (|201ll ') presented an anal- 
ysis of a z = 1.82 QSO associated with submm 
emission and which is surrounded by an over- 
density of submm galaxies (both QSOs and submm 
galaxi es are strong tr ac ers o f large scale st ru cture , 
e.g. ISmith fc Heckmaiil Il990l: lEUirigson et all Il99ll: 



Austermann et al. 20091: IChapman et al. 20091: Bielbv et al] 



2OIOI : ISiemiginowska et all boid H Matsuda et al. I l201lh . 

Based on spectral energy distribution (SED) fltting, they 
conclude that the emission in the mid-infrared (MIR) 
range is dominated by reprocessed AGN emission, whilst 
the submm is dominated by a starburst co ntribution. A 
simila r result was found for a z ~ 4 QS O bv lCoppin et al.l 
l|2009t ). However, iHill fc Shank3 (|2011al ) fit galaxy counts 
and colours using optically deflned Pure Luminosity 
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Evolution (PLE) models, where dust reradiates absorbed 
optical light into infrared spectra composed of local galaxy 
templates, across the MIR bands from 3.6/im to 8/im and 
up to a redshift of at least z = 2.5. They show that a 
significant contribution from AGN to the 250 — 870 /xm 
SED would remove the need to invoke a top-heavy IMF for 
high-redshift starburst galaxies. It may therefore be that 
reprocessed emission from AGN could be more important 
than previously thought. 

Interestingly, a significant AGN contribution to the 
submm population could have relevance for the origin of 
the hard X-ray (i.e. 10-30 keV) background. At softer en- 
ergies (i.e. < 10 keV), much of the X-ra y background 
(XRB ) is resolved into sources, usually AGN (|Shanks et al.l 
Il99lh . At harder energies therefore the expectation is that 
heavily obscured AGN may form t he background (e.g. 
IComastri et al.lll995l: iGilh et all 120071 ). However, the hard 
X-ray sources that have so far been identified by Chan- 
dra and XMM are skewed towards lower redshift, peak- 
ing at 2 < 1 and hence have so mewhat lower luminos i- 
ties, Lx « lO''^ - 10** erg/s (e.g [Alexander et all [200^ . 
Most studies conclude that the contributing sources so far 
detected only comprise < 25% of the XRB (IWorslev et al ' 



20051; iKrivonos et~aLll2005l: IWor slev et "aLll2006l : iDaddi et al 



20071 : iHickox fc Markevi tch 20o3). Should submm sources 



contain obscured AGN, this may go some way to explaining 
the missing XRB flux. 

There is some evidence that different populations of 
high redshif t galax ies host obscured AGN. For example, 
IPaddi et al.l ([2003) reported that a number of BzK se- 
lected galaxies show a thermal, warm dust excess at 24 
microns resulting from faint, hard X-ray sources. The sug- 
gestion is that these are obscured QSOs, but again these 
X-ray sources may be too faint to account for any more 
than 25% of the backgroun d at 10-30 keV (whilst we note 
that [Alexander et al.l I2OIII have recently reanalyzed the 
iDaddi et all |2007| galaxy sample with deeper X-ray data 
and found a factor of > 10 fewer obsc ured AGN than orig- 
inally suggested bv lDaddi et al.l [20071 ). In terms of submm 
galaxies, stacked X-ray spectra show a broad Fe Kc, line 
l| Alexander et al1l2005al ). whilst many of the optical spec- 
tra of submm galaxies show evidence for broad lines. Ad- 
ditionally, high ionisation lines indicative of AGN activ- 
ity have also been detected at FIR wavelengths in Her- 
schel SPIRE FTS obse rvations of a z ~ 3 submm galaxy 
IValtchanov et al]|201ll ). Finallv. iBusswell fc Shanksl (I2OOII ) 
suggested that the obscured QSOs might explain the bright, 
«5mJy, submm number counts if the QSOs had a temper- 
ature of Ri30K (see also lAlmaini et al] 1 19991 ) and this is at 
least consistent with the temperatures being reported for 
the SCUBA Half- Degree Extrag a lactic Survey (SHADES) 
submm sources bv ICoppin et al. (|2008bh . Thus, it may be 
that submm galaxies account for a further fraction of the 
X-ray background at 10-30 keV. 

We have undertaken a survey designed to measure 
870^m fluxes for a sample of known quasars with the aim 
of comparing the properties of obscured versus unobscured 
AGN. This has be en performed in the William H erschel 
Deep Field rWHDF: [Metcalfe et al ]|l995l . l200ll . l2006t ). which 
is especially suitable for this survey thanks to a signifi- 
cant sample of spectroscopically confirmed quasars within 
a small, easily observable area. This includes both unob- 



scured and heavily obscured sourc es matched in redshift and 
luminosity (| Vallbe-Mumbrull2004l ) . Although not capable of 
producing definitive results on their own, the observations 
we present here are a crucial step towards identifying the 
contribution of AGN to the production of the submm emis- 
sion in submm sources and thus, this population's contri- 
bution to the submm and X-ray backgrounds. This paper 
therefore presents the results from the submm survey, fo- 
cu ssing on statist i cal ana lyses that complement the analyses 
in lHill fc Shank j l|2011bl ). and which will form the basis for 
further multi-wavelength treatment of this issue. 



2 THE WILLIAM HERSCHEL DEEP FIELD 

The William Herschel Deep Field is a 7'x7' area centred 
at 00''20'" -1-00° ( J2000) which has a wealth of multiwave- 
length dat a and has been extensively studied over the pa st 
15 vears (iMetcalfe et all Il995l : iMcCracken et all l2000al lbl: 
iMetcalfe et al.ll200ll , 12009 ). 

The field has ultradeep, ground-based optical UBRIZ 
imaging from the William Herschel Telescope in La Palma 
(reaching B < 27.9) as well as near-infrared H and K imag- 
ing from Calar Alto and the UK Infrared Telescope and 
deep, high-resolution J-band imaging from the Hubble Space 
Telescope's Advanced Camera for Surveys. 

In addition to this comprehensive optical/NIR cover- 
age, the WHDF also has Chandra X-ray coverage, reach- 
ing a depth of « 10~^^ ergs s~^ cm~'^ over the whole area 
with a total integration time of ~ 70 ksec. These obser- 
vations were undertaken between November 2000 and Jan- 
uary 2001. 170 X-ray sources were detected at ^ 2a signif- 
icance, of which 69 were at ^ 3a and 36 at ^ 5a. Spec- 
troscopic follow-up of some of these X-ray sources was sub- 
se quently performed and these observations are described 
bv lVallbe-Mumbrul (|2004l ) and in t]5. II below. Further spec- 
troscopic data (of star- forming and passive galaxies r espec - 
tive ly) has also been presented by iBohm fc Zieglerl (|2007n 
and lFritz et"al] (|2009t ). 

The most recent additions to the WHDF data are (o) 
an 870/im submm survey, which we undertook between 2008 
and 2009 and which is described in this paper, and (6) radio 
observations at 8.4 GHz (A — 3.6cm) acquired with the Ex- 
panded Very Large Array (EVLA) in New Mexico in 2010. 
The reduction of these radio data is under way. 



3 SUBMM OBSERVATIONS AND DATA 
REDUCTION 

3.1 Observations 

We have acquired 21 hours of observations of the WHDF 
with the Large A pex Bolometer Camera (LABOCA; 
ISiringo et al] 120091 ) on the 12 m APEX telescope 
(|Giisten et al.|[2006l ). The LABOCA instrument comprises 
295 semiconducting composite bolometers arranged in a se- 
ries of concentric hexagons. LABOCA is sensitive to radi- 
ation in a passband centred at 870/im, with a FWHM of 
« 150/im. 

The LABOCA beam has a FWHM of 18. 6" and the total 
field of view (FoV) of the detector is 11.4'. The WHDF cov- 
ers a region of « 7'x7', as noted above, so this field fits well 
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R.A. (J2000) 

Figure 1. (a) LABOCA 870^m intensity map of the WHDF, with 11 sources detected at ^ 3.2a- marked by open blue circles. The 
overlayed contours show the redcued noise map, with contour levels marked in mjy/beam. Spectroscopically identified QSOs (listed in 
Tabled are also plotted, coded by their estimated hydrogen column density. The reduced ofF-axis sensitivity of LABOCA produces the 
increase in noise at the edges of the frame, (b) LABOCA signal-to-noise map of the WHDF, with the same 11 sources marked. Again 
the contours show the reduced noise map. 



into the FoV of LABOCA. The LABOCA detectors do not 
form a contiguous array, however, so to achieve full sampling 
of the field we carried out our observations using a standard 
spiral raster map pattern. In this process, the centre of the 
array is moved in a spiral pattern and whilst being shifted 
laterally in a raster configuration, in order to fully sample 
the FoV. 

Our observations were carried out in two separate ob- 
serving runfl, the first on 29-30th August 2008 and the sec- 
ond on ll-12th May 2009. The observing conditions were 
very good, with the range of zenith opacities, r^, spanning 
0.17 to 0.23 in the first run and 0.17 to 0.28 in the sec- 
ond (these being calculated as standard from a linear com- 
bination of LABOCA 870 fiui skydip results and opacities 
determined from the APEX radiometer). 

3.2 Data reduction 

Initial data reducti on was performed usin g the standard BOA 
pipeline software (jSchuUer et al.l l201Cf I) using the recom- 
mended sequence for 'weak sources'. Counts from the de- 
tector are first converted into an output voltage, which can 
then be converted into a flux density, in standard units of 
Jy beam~^, via a voltage-flux relation determined empir- 
ically during LABOCA's commissioning. Sky removal was 
performed in BoA using the iterative medianNoiseRemoval 
command, which computes and corrects for the relative 
gains for a channel with respect to the mean signal. In addi- 
tion to the usual de-spiking, data flagging, correlated noise 
removal and Fourier space filtering during the reduction, a 
DC offset (modelled by a first order polynomial) is removed 

1 ESO programme IDs 081.A-0897(A) and 083.A-0707(A) 

2 www. mpifr-bonn. mpg.de /div/ submmtech / software /boa /boaman. 



from each scan for each bolometer. A further correction is 
then applied to account for the zenith opacity at the time of 
observing. The zenith opacity, r^, is a, measure of how much 
incoming radiation is absorbed by the atmosphere and is de- 
termined by 'skydip' observations (calibration exposures of 
the sky). Flux calibration was performed using observations 
of primary and secondary calibrators, with Uranus and Nep- 
tune being used as primary calibration sources during our 
observing run. The result of this process is a series of expo- 
sure maps for each scan of the target field. The maps were 
then co-added using BoA, which weights the signal from 
each scan image by l/a^, where a is the pixel-noise and is 
estimated by adding in quadrature the noise levels for each 
bolometer that contributes to a given pixel. The final output 
is a combined intensity (fiux) map, as well as corresponding 
rms noise and signal-to-noise (SNR) maps. Each of the maps 
has a pixel scale of 9".l/pixel. 

This data reduction process was carried out separately 
for the data from each of the two observing runs. As a final 
stage in the reduction process, therefore, we combined the 
maps, weighting them by the noise. The final intensity map 
was produced according to Eqn. [T] and the final noise map 
according to Eqn. (2] 
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where I is intensity, a is rms noise and subscripts 1 and 
2 indicate the first and second observing runs. A signal- 
to-noise map was produced by taking the ratio of the two, 
and each of the maps was then Gaussian-smoothed using 
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the 18". 6 beam profile, which gave a final resolution for the 
smoothed m aps of 27" . The sm oothing was comparable to 
that done bv lWeifi et ai] (|2009l '). who have shown that flux 
estimates from similarly smoothed maps are consistent with 
non-smoothed fluxes. 



4 LABOCA DATA 

The final intensity and SNR maps are shown in Fig. [T] It 
is clear from Fig. [T] that the noise increases substantially at 
large off-axis angles; in the SNR map (right-hand panel), this 
increase in the noise manifests itself as a dearth of sources 
around the edges of the frame. In the central 100 arcmin^ 
area, the maps are very similar, with points of bright inten- 
sity having corresponding peaks in the SNR map, an indica- 
tion that the noise level is relatively uniform across the field 
centre. 

4.1 Noise level 

Across the central 100 arcmin^ region of the field the data 
reach an rms noise level of < 2 mjy beam~^, making the 
WHDF one o f the deepest sub mm fields observed to date 
(c/. figure 6 of IWeifi et aLll2009l ). 

To check the reliability of the pipeline-reduced noise 
map we have produced a 'standard deviation map', by mea- 
suring the standard deviation of the intensity map in a series 
of annuli from the field centre, having masked out the bright 
sources (see Section 14. 2|) . If the pipeline has worked success- 
fully, this standard deviation map should be comparable to 
the BoA-produced noise map. 

In Figs. [2] and |3l we compare the contours and radial 
profiles, respectively, of the two maps. For the pipeline- 
reduced map, the noise profile shown in Fig. [3] is the me- 
dian of 120 profiles measured radially at 3' intervals. The 
agreement between the two noise maps is very good — this 
is particularly clear from Fig. [3] — suggesting that the noise 
has been reliably estimated for our field and that our quoted 
depth of 2 mJy beam~^ over an area of 100 arcmin^ is ro- 
bust. 

4.2 LABOCA sources 

The field contains 11 significant submm sources, circled in 
Fig.[T] We select the sources from the SNR map, with a cri- 
terion of S/N ^ 3.2. This significance was chosen by com- 
paring the map to its inverse: there are no negative spikes 
in the SNR map with a magnitude of 3.2it. This is illus- 
trated in Fig. 21 which shows a pixel value histogram for 
the SNR map, with a Gaussian fit (blue curve). The map 
shows a strong positive excess while the negative side of the 
distribution follows the Gaussian curve well. Based on the 
Gaussian fit, we estimate that the data would contain « 0.2 
spurious peaks above our chosen SNR limit within the cen- 
tral 100 arcmin'^ region. 

The positions, SNRs and fiuxes of the 11 sources are 
summarised in Table [T] We note that the quot ed fluxes have 
not been corrected for the flu x bias effect (|Coppin et al.l 
l2008al : lAustermann et al]|2009l l , however we do not consider 
this necessary for the scope of this paper. The positions of 
the 11 sources are indicated on Figs. [T] and [2] and their fiuxes 
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Figure 2. Noise contour maps of the WHDF, with the rms noise 
level marked at 1.2, 1.5, 2.0 and 5.0 mJy beam~^, using (a) the 
pipeUne-reduced noise map and (h ) the standard deviation of the 
background flux measured in radial annuli. The image pixel size 
is 9."1. 3.2(T LABOCA sources are marked as blue circles. The 
maps generally agree well (see also Fig.[3ll. Most of our LABOCA 
sources are detected within the central region where the noise is 
lowest. 



T 




1 r / 

Pipeline noise map 

Stdev noise map 



1 1 1 I 1 1 1 I 1 1 1 I 1 1 L 

2 4 6 

r (arcmin) 

Figure 3. Radial noise profile of the WHDF showing the rms 
noise level, as a function of radius from the centre, for both the 
pipeline-reduced and standard deviation noise maps, described in 
the main text. 3.2(t LABOCA sources arc marked (blue circles) 
at their respective source flux and radial distance from the centre 
of the field. 

are shown in comparison to the background noise in Fig. |3l 
Most of the sources lie within the < 2 mjy/beam central 
area of the field, with 7/11 located within the a ^ 1.2 mJy 
contour. 



4.3 Completeness 

We estimate the completeness of our submm observations 
using simulated sources placed in the final submm map 
across a range of fluxes. At each fiux level, 500 sources are 
individually added to the map in turn and the detection 
analysis re-performed for each source. This is limited to the 
< 2 mJy/beam region only and excluding regions within 2 
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Table 1. WHDF LABOCA 870 fim sources detected at 3.2a. 



ID 


RA 


Dec 


SNR 


5870 °- 


A0 




(J2000) 




(mjy/beam) 


(") 


WHDF-LAB-01 


00:22:37.55 


+00:19:16.8 


4.3 


5.1 


9.4 


WHDF T AR n? 

V V -L J- l-J J- -Li ii J_j U 


00:22:28.44 


+00:21:42.6 


4.2 


4.3 


9.6 


WHDF-LAB-03 


00:22:46.06 


+00:18:40.3 


3.8 


5.4 


10.7 


WHDF-LAB-04 


00:22:29.66 


+00:20:20.5 


3.7 


4.1 


10.9 


WHDF-LAB-05 


00:22:22.97 


+00:20:11.4 


3.6 


4.0 


11.3 


WHDF-LAB-06 


00:22:32.09 


+00:21:24.3 


3.6 


3.9 


11.3 


WHDF-LAB-07 


00:22:48.49 


+00:16:32.8 


3.5 


8.2 


11.6 


WHDF-LAB-08 


00:22:29.66 


+00:16:05.4 


3.4 


6.2 


11.9 


WHDF-LAB-09 


00:22:19.90 


+00:17:00.1 


3.2 


5.1 


12.7 


WHDF-LAB-10 


00:22:35.16 


+00:24:08.3 


3.2 


3.3 


12.7 


WHDF-LAB-11 


00:22:25.40 


+00:20:11.4 


3.2 


3.4 


12.7 



" we note that listed fluxes are measured fluxes and not corrected for flux bias. 
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Figure 4. A histogram of pixel values from the SNR, map. A 
Gaussian has been empirically fit to the profile — most of the 
pixels lie within this distribution, however a significant positive 
excess appears due to the presence of sources. The shaded region 
indicates those pixels lying above the 3.2(t threshold we have set 
for source detection. 

pixels of any of the 11 detected sources. The result is shown 
by the filled black circles in Fig. [S] with the curve showing 
a polynomial fit to the result. Based on this analysis, we 
find the data to be 50% complete at the 4.3 mJy/beam level 
(shown by the dotted line) . 

5 SUBMM OBSERVATIONS OF WHDF 
QUASARS 

5.1 Obscured and unobscured quasar samples 

In 321 we briefly summarised the Chandra observations of 
the WHDF, which yielded detections of 170 X-ray sources. 




0.0 I , , , I , , , I , , , I , , , I , , , 

2 4 6 8 10 

(mJy/beam) 



Figure 5. The completeness of the submm data (< 2 mJy/beam 
region only) estimated based on simulated sources placed in 
the submm image. The solid line shows a polynomial fit to the 
completeness estimates as a function of source brightness, with 
50% completeness marked by the dotted line at Sgro/jm = 4.2 
mJy/beam. 

In 2001 and 2002, 36 of these WHDF X-ray sources were 
targetted in 10 hours of spectroscopic follow-up with the 
LDSS2 multi-object spectrograph on the 6.5m Magellan-1 
(Walter Baade) telescope at the Las Campanas Observatory 
in Chile. These were optical observations, using a grism cen- 
tred at A ^ 5500A. 

A stated aim of the LDSS2 observations was to in- 
vestigate non-quasar X-ray sources, so the principal tar- 
gets for spectroscopic follow-up were X-ray sources whose 
optical counterparts did not appear especially point-like 
(j Vallbe-Mumbrull2004l ) . Nevertheless, 15 of the targets were 
spectroscopically conflrmed as QSOs, and a further 2 were 
classified as being either QSOs or narrow emission line galax- 
ies (NELGs). In this section we take the 15 confirmed QSOs 
as our sample of WHDF quasars; details of these sources are 
given in Tabled 

The QSO classification was made by IVallbe-Mumbrd 
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Table 2. Spectroscopically confirmed quasars in the WHDF. The positions, 0.5-10 keV fluxes (in ergs s ^ cm X-ray hardness ratios 
and spectroscopic redshifts are given, obscured quasars are indicated in bold. 



ID" R.A. Dec. 5o.5-io HR z 

(J2000) (erg/cm2/s) 



WHDFCH005 


00:22:35.963 


00:18:50.04 


5.62 


X 


10 


-14 


—0.60 


± 


0.08 


0.52 


WHDFCH007 


00:22:24.821 


00:20:10.94 


1.17 


X 


10- 


-14 


-fO.82 


± 


0.33 


1.33 


WHDFCH008 


00:22:22.884 


00:20:13.24 


3.62 


X 


10- 


-15 


-0.20 


± 


0.24 


2.12 


WHDFCH016 


00:22:45.164 


00:18:22.64 


1.44 


X 


10- 


-14 


-0.43 


± 


0.14 


1.73 


WHDFCH017 


00:22:44.468 


00:18:25.64 


3.22 


X 


10- 


-13 


-0.55 


± 


0.03 


0.40 


WHDFCH020 


00:22:36.142 


00:24:33.84 


1.09 


X 


10- 


-14 


-0.55 


± 


0.17 


0.95 


WHDFCH036 


00:22:31.734 


00:25:38.84 


6.26 


X 


10- 


-14 


-0.48 


± 


0.06 


0.83 


WHDFCH044 


00:22:55.092 


00:20:55.74 


2.66 


X 


10- 


-14 


+0.60 


± 


0.14 


0.79 


WHDFCH048 


00:22:41.297 


00:25:33.34 


2.15 


X 


10- 


-14 


-0.43 


± 


0.10 


1.52 


WHDFCH055 


00:22:11.862 


00:19:50.44 


2.17 


X 


10- 


-14 


-0.23 


± 


0.06 


0.74 


WHDFCH090 


00:22:48.795 


00:15:18.74 


4.83 


X 


10- 


-14 


-0.41 


± 


0.06 


1.32 


WHDFCH099 


00:22:11.187 


00:24:04.13 


8.84 


X 


10- 


-15 


+0.11 


± 


0.18 


0.82 


WHDFCH109 


00:22:09.917 


00:16:28.94 


6.69 


X 


10- 


-14 


-0.55 


± 
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Figure 6. Optical spectra of WHDFCH007 (top) and WHDFCH008 
(bottom) acquired with LDSS2 on the Magellan telescope. Both 
objects are associated with submm sources in our sample. Broad 
Mgll emission is observed in the spectrum of WHDFCH007. Only 
narrow emission lines are d etected for whd f chOO S, which was 
classified as a type 2 QSO bj(wibi^Mumb3 l|2004l l. 



(|2004l ) on the basis of both X-ray luminosity (Lx > 10*"* ergs 
s~^) and optical emission lines. Most of the sources classed 
as QSOs showed broad emission lines, including Nev A1240, 
Siiv A1400, Niv A1486, Heii A1640, Oiii] A1663 and/or Niii 
A1750, all of which are AGN indicators. If broad lines were 
detected the source was classified as a type 1 quasar, for 
example WHDFCH007 which is shown in the top panel of Fig. 
El One source was classified as a type 2 quasar on the basis 
of strongly detected narrow emission lines (whdfch008); its 
spectrum is shown in the bottom panel of Fig. |6l 

The level of obscuration of an AGN can be estimated 
directly from the X-ray spectrum using the hardness ratio. 



HR. This is defined as HR = {H~S)/{H + S), where H and 
S represent the photon counts in the hard (2-8 keV) and 
soft (0.5 - 2 keV) X-ray bands respectively. Increasing lev- 
els of absorbing hydrogen column density produce greater 
absorption in the soft band than in the hard, so leading 
to larger values of HR. Assuming a given intrinsic spectral 
slope, F, the absorbing hydrogen column density can be es- 
timated. Fig. [7| shows hardness ratio against redshift for our 
quasar sample, compared to the predicted tracks for quasars 
at different absorbing column densities (an intrinsic photon 
index of F = 2 is assumed). 

The locus of the quasars is at HR ~ —0.5, consistent 
with the model that predicts HR = —0.5 for essentially all 
unobscured QSOs. Only four of the 15 QSOs are harder than 
HR = —0.2; these 4 sources — WHDFCH007, -008, -044 and 
-099 — are expected to be highly obscured. whdfchOO? and 
WHDFCH044 have extremely hard spectra with HR > 0.6, 
corresponding to an apparent photon index of F ^ —1. Such 
sources are relatively rare, with for example only 2 HR « 
0.6 objects being reported in the ~ 2 deg^ XMM-Newton 
COSMOS survey jMainieri et all bOOTl ) (which we note is 
less deep by a factor of « 2 — 3x than the WHDF X-ray 
data used here, whilst they do not present HR values with 
errors of > 0.3). 

Based on Fig. [7| we divide the WHDF quasar sample 
into three groups: heavily obscured (A^h > lO'^^ cm^^), 
mildly obscured (A^h ~ 10^^ cm^^) and unobscured (A^h < 
10^^ cm"'^); these are marked in both Fig. [1] and Fig. [7|by 
black, blue and cyan stars respectively. 



5.2 submm properties 

5.2.1 Possible counterparts 

We now look for associations between the submm sources 
found in the LABOCA data and the spectroscopically con- 
firmed QSOs described above. For the purposes of this anal- 
ysis, we restrict the sample to only those QSOs that lie 
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Figure 7. Hardness ratio versus redshift. The curves show pre- 
dicted tracks for obscured QSOs at different column densities, 
which are indicated on the figure (in units of cm~^). An intrin- 
sic r = 2 power-law spectrum is assumed. Below TVjj = 10^^ 
cm~^ the lines become indistinguishable. WHDF X-ray QSOs 
are marked; each is labelled with its ID (with the WHDFCH prefix 
omitted). On the basis of this figure we classify the QSOs into 
3 categories: heavily obscured (A'^jj > 10^^ cm~^; black), mildly 
obscured (Ah ~ 10^^ cm^^; blue) and unobscurcd (Ah < 10^^ 
cm~2; cyan). 

within the 2 mjy/beam contours of the pipeline reduced 
submm noise map shown in the left hand panel of Fig. [2l 
This leaves two Nh > 10^^ cm"^ three Nh ~ 10^^ cm~^ 
and four Nh < 10^^ cm"^ QSOs. 

In order to identify coincident sources, we use a max- 
imum separation between a QSO and a submm source of 
2.5A9, where A6 is the positional uncertainty on the given 
submm source (listed in Table [T]). Given the signal-to- noise 
hmit on our submm catalogue of S/N — 3.2a, this gives a 
maximum po ssible separation of 2.5(0.66I(5'/A')"^) = 12.7" 
(|lvison et al.ll2007l ). with the smoothed map resolution of 
6 = 27". The majority of the submm sources will have a 
separation constraint somewhat smaller than this however, 
given their higher signal-to-noise. 

In order to evaluate the significance of a given align- 
ment, we use the corrected Po isson probability, P, as em- 
ployed bv iDownes et all (| 19861 '). As such, we estimate the 
probability of a chance alignment using the observed inte- 
grated sky density of X-ray sources, N{> S), as a function 
of the soft X-ray flux, 50.5 -2, which we calculate based on 
the power-law fit given bv IVallbe-Mumbrul (|2004) for the 
WHDF sources: 

logAf(> So.5-2) = -8.6 - 0.761og(So.5-2) (3) 

where A'^(> So. 5-2) is in units of deg^^ and So. 5-2 is in 
units of ergs/cm^/s. We note that this is in good agree- 
ment with other such f its to the X-ray num ber counts (e.g. 
iMushotzkv et aL l l200(]| : lGiacconi et alj|200ll ). 

Based on the chosen 2.5A9 limits for identifying coinci- 
dent sources, we find that two QSOs have possible LABOCA 
counterparts: WHDFCH007 (coincident with whdf-lab-11) 
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Figure 8. Thumbnail images of the LABOCA SlOfim intensity 
map at the positions of the highly obscured {Nfi > IQ-^-^ cm^^) 
QSOs WHDFCH007 and whdfchOOS. In each case the X-ray source 
position is marked by the star and nearby submm sources are 
marked by the blue circles (with the radii giving the estimated 
2.5A9 positional accuracy). Both sources have closely associated 
submm emission. 



and WHDFCHOOS (coincident with whdf-lad-05). Fig. [8] 
shows thumbnail images of the LABOCA flux map at the 
positions of these X-ray sources, and > 3.2a submm sources 
are seen close to the two QSOs. In each case, the star 
shows the location of the QSO and the blue circles show the 
2.5A6' limit around nearby submm sources. The corrected 
Poisson probabilities that these are chance alignments are 
P = 0.029 and P = 0.004 (i.e. 2.9% and 0.4%) for whd- 
FCH007 a nd whdfchOOS res pectively. For comparison, we 
note that llvison et all (|2007l ). use a limit of P ^ 0.05 as a 
constraint for secure alignments between submm and 24/im 
sources. 

These two alignments are with QSOs in our highly ob- 
scured sample, whilst we note that the other two highly 
obscured quasars, WHDFCH044 and -099, lie outside the im- 
posed 2 mjy beam~^ noise limit. It is interesting, therefore, 
that of the two highly obscured WHDF QSOs which lie in 
the central, low-noise area of the field, both appear to have 
bright submm counterparts. 

In Figs. [9] and IIOI we show thumbnail images for three 
Nn ~ 10^^ cm~^ QSOs and four unobscured QSOs, respec- 
tively. Again we note that the full samples have 6 and 8 
sources respectively, but those QSOs not shown are rejected 
due to lying outside the < 2 mJy beam~^ noise region. 

Some of the sources in these figures appear to lie close 
to bright areas, e.g. WHDFCH016 in Fig. [9]or WHDFCH017 in 
Fig. 1101 however, none of these QSOs could be said to be 
coincident with a peak in the map as was the case for whd- 
FCH007 and -008, lying as they do well outside the beam po- 
sitional accuracy limits. We note that at the fluxes of the two 
submm sources associated with obscured AGN (i.e. ~ 3 — 4 
mJy/beam), our data is ~ 30% complete, which could com- 
promise the observation that the Nh ^ lO^'^ cm~^ QSOs 
are not associated with submm emission to the level of the 
Nh ^ 10^^ cm~^ QSOs. In the following section we there- 
fore perform a stacking of each of the AGN populations to 
look for signatures of submm sources close to our detection 
threshold. 
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Figure 9. As in Fig. |8] but for sources characterised as mildly 
obscured with Afjj 
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Figure 11. Noise-weighted stacked intensity maps for the three 
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Figure 10. As in Fig. |8] but for sources characterised as unob- 
scured with A^h < 10^^ cm~^. 



5.2.2 Stacking 

To get an overall picture of the submm flux associated with 
our highly obscured, mildly obscured and unobscured quasar 
populations, we stack the submm flux maps for each of the 
populations. To ensure that the noisy fringe sources do not 
dominate the stacked flux there are two ways to proceed 
— either to exclude the sources near the edge and stack 
only the sources linearly within the < 2 mjy/beam region 
(exploiting the relatively uniform noise across the centre), or 



imum search radius of = 12". 7 used in matching the QSOs 
and submm sources. Only the most highly obscured QSOs show 
significant stacked submm emission. 



to perform a noise- weighted stack which mitigates the effect 
of the noisy objects. Having performed both analyses, we 
find the results are not significantly affected by the choice 
of method. In the analysis presented here, for simplicity we 
use stacks incorporating only those sources within the < 2 
mJy/beam region, without any noise weighting. 

The stacked flux images for the three QSO popula- 
tions are shown in Fig. 1111 incorporating two sources in the 



Nh > 

-2 



10 cm sample, three sources in the Nh 



10^ 



cm sample and four sources in the Nh < 10 cm^ sam- 
ple. Taking the peak value within our maximum search ra- 
dius of 12. "7, we find a peak signal of Ssroum = 3.4 ± 0.8 
mJy/beam at an angular distance of 5" from the QSO po- 
sitions for the Nh > 10^^ cm"^ stack. For the Nh « 10^^ 
cm~^ and Nh < 10^^ cm^^ stacks we find peak signals 
of S'sTOMm = 1.9 ± 0.8 mJy/beam and Sstom™ = 0.7 ± 0.5 
mJy/beam respectively (both being found at an angular dis- 
tance of ~ 10" from the stack centre). These plots reinforce 
the observations based on individual X-ray sources, that the 
highly obscured quasars are the only sources to be associated 
with submm bright sources, without being affected by the 
completeness of the submm data at the ~ 3 — 4 mJy/beam 
flux level. 



5.2.3 Statistical power 

It is important to account for the statistical power afforded 
by the small numbers of objects involved in the current 
analysis. We therefore estimate the confidence limits of the 
sample size presented here (based o n the tabulate d small 
sample Poisson statistics provided bv lGehrei^ll986l ) for the 
fraction of unobscured and obscured AGN associated with 
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submm emission. Based on the sample of 2 obscured AGN, 
we estimate that our observations give a lower limit of 35% 
on the percentage of such objects having associated > 3 
mjy/beam submm emission (at the la level). Conversely, 
the sample of Nh < 10^^ cm~^, gives an upper limit of 37% 
on the percentage of these unobscured AGN being associ- 
ated with > 3 mJy/beam submm emission (again at the la 
level). Although the samples considered here are small and 
cannot fully constrain the relationship between AGN obscu- 
ration and submm emission, our analysis remains consistent 
with a model in which submm emission arises preferentially 
from more obscured AGN - in agreem e nt with the analy - 
ses of, for exa r nple. iPage et all (|2004l '): iLutz et all (|2010f ): 
iHiU fc ShanQ (|2011bh . 



5.3 Future improvements 

Despite the relatively poor resolution of our submm data, it 
appears that each of the X-ray QSOs WHDFCH007 and whd- 
FCH008 has only a small chance of accidentally aligning with 
a submm source. It is notable that only the most obscured 
AGN show possible 870^m detections and these include the 
Type 2 QSO, whdfchOOS. Stacking the submm data as a 
function of X-ray absorption also provides further evidence 
that faint, X-ray obscured QSOs are preferentially submm 
bright. 

Further tests of the reality of the associations between 
the two X-ray obscured QSOs and the LABOCA sources will 
soon be available from the 8.4GHz EVLA radio survey of the 
WHDF, which will have ~ 2x improved spatial resolution 
compared to the LABOCA data. These data are still at the 
reduction stage and the results will be reported in a later 
papei0. 

In addition, we are seeking Herschel observations at 
250, 350 and 500/im, which will constrain the SED and allow 
dust temperatures to be measured for the detected sources. 
If the unobscured QSOs are found to be detected in these 
shorter- wavelength bands, it could reveal whether any tem- 
perature difference exists betwee n obscured and unobscured 
quasars, as was hypothesised bv lHill fc Shanks! (|2011bl ). 

Finally, we are also proposing to observe the two 
X-ray obscured QSOs with Atacama Large Millime- 
tre/submillimetre Array (ALMA) in its extended configura- 
tion to investigate further the submm counterparts of these 
sources. Not only will ALMA give better positional informa- 
tion but its « 0".5 angular resolution will allow us to test if 
the associated submm sources are point-like or extended at 
< 3kpc spatial resolution. 



6 CONCLUSIONS 

We have presented an 870/im survey of the well-studied 
William Herschel Deep Field, reaching a depth of < 2 
mJy/beam over an area of 100 arcmin'^. In total, 11 sources 
have been detected at a significance of ^ 3.2cr. From the 
noise distribution of the image, we estimate the number of 
false detections in the 100 arcmin^ region to be ~ 0.2, whilst 



^ Wc note that since submission of this paper, both WHDFCH007 
and whdfchOOS have been detected in the 8.4GHz EVLA data. 



the completeness in the same region is 50% at a flux of 4.3 
mJy/beam. 

We find that 2 of the 11 submm sources are likely coun- 
terparts of X-ray-selected WHDF quasars: WHDFCH007, an 
extremely hard QSO with F fa — 1 and whdfchOOS, a type 
2 QSO with only narrow optical emission lines, both con- 
sistent with being heavily obscured. We divide a sample of 
15 WHDF quasars into three subsets, having A^h < 10^^, 
A^H ~ 10^^ and A^n > 10^^ cm"^, and find that only 
the most obscured population shows any significant stacked 
870Atm flux. 

Although based on a relatively small sample, our find- 
ings are supportive of a model in which obscured AGN are 
submm bright and unobscured AGN are not. This picture 
of faint X-ray obscured AGN being preferentially stronger 
submm emitters is difficult to fit into the unified AGN 
model without some adaptation, given that both X-ray 
obscured and unobscured AGN would be expected to be 
equally strong submm emitters in this case. Never t heless , 
the observationa l results of, for example. | Page et al] (|2004[): 
lAlexarider et al.1 (|2005al ): iLutz et al.1 (|2010l ): iHill fc Shankd 
(|2011bl ) and now the results presented here in the WHDF 
all tend to support the idea that obscured AGN are more 
frequently submm emitters than unobscured AGN. 

From these results, two solutions have been proposed: 
either (a) AGN heating is a dominant component of dust 
heating, producing submm emission (e.g. Granato fc Danese 



1994 



Willott et al]|2002l : iGrimes et al.ll2005l : iHill fc Shank: 



2011bj) or (b) obscured AGN are an evolutionary concurrent 
phenomenon with strong star-formation and it is the star- 
formation that is the dominant mechanism for the dust heat- 



ing in these objects (e.g. Page et al. [2004 : 1 Alexander et al.l 



l2005bl : Ivignah et al.H2009l: iLutz et al.ll2010ir 

Assuming the former. llliu fc Shanksl (|2011bh argue that 
the contribution of obscured AGN to the submm background 
may be as much as f» 40% (compared to « 13% in the uni- 
fied/evolutionary model case). In addition, the same authors 
suggest that, should the primary driver of submm emission 
in these sources be heating from obscured AGN, then these 
may be the source of a significant fraction of the > 10 keV 
XRB, whilst the contribution would be much lower should 
heating from star-formation be the dominant mechanism. 
Other potential candi dates still fall short of fully account- 
ing for the XRB (e .g. IWorslev et a"l]|2005l : iGilli et al.|[2007l : 
iTreister et al.|[2009l ) and in this context, quantifying the con- 
tribution from obscured AGN in the submm source popula- 
tion is an important goal. Although we are unable to inform 
fully on this based on the results here, these observations 
remain an important step in the process. 

Upcoming radio data, as well as proposed observations 
in the submm and other bands, will enable further investi- 
gation of the submm properties of obscured and unobscured 
AGN in the WHDF. As discussed, Herschel observations 
will add important constraints to the SEDs of these ob- 
jects and facilitate the measurement of dust temperatures, 
whilst the obscured QSOs will also make excellent targets 
for high resolution observations with ALMA. Crucially, such 
observations would not only improve on the accuracy of the 
submm source positions, but would also inform on whether 
the submm sources are extended, indicating a galaxy- wide 
starburst origin, or point-like and therefore more associated 
with the activity in the nucleus. 
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